TRAJECTORY OPTIMIZATION AND GUIDANCE 
FOR AN AEROSPACE PLANE 


l/j/Jic e/ 


Final Report 

September 1988 through August 1989 


/jj- 05' 

J3 / 7-2 <3 


September 1989 



Research Supported by NASA Langley Research Center 
NASA Grant Number NAGI-907 


Principal Investigator: Kenneth D. Mease 
Graduate Research Assistant: Mark A. Van Buren 


Department of Mechanical and Aerospace Engineering 
School of Engineering and Applied Science 
Princeton University, Princeton, New Jersey 08544 


NASA Grant Monitor: D. D. Moerder 


(NASA-CR-l 85884 ) TRAJECTORY OPTIMIZATION 
AND GUIDANCE FOR AN AEROSPACE PLANE Final 
Report, Sept, i960 - Aug. 1989 ( p rinceton 
Univ.) A2 p CSCL OIC 

G3/05 


N90-13369 


Unci as 
0231728 


TABLE OF CONTENTS 


ACKNOWLEDGEMENTS 2 

NOTATION 3 

SUMMARY OF WORK COMPLETED 4 

1. INTRODUCTION 5 

2. PROBLEM FORMULATION 6 

2.1 Equations of Motion 6 

2.2 Initial and Final Conditions 7 

2.3 Performance Index 8 

2.4 State and Control Constraints 8 

2.5 Optimality Conditions 8 

3. VEHICLE MODELING 12 

3.1 The Aerodynamic Characteristics 12 

3.2 The Propulsion System 12 

3.1.1 The Ramjet 13 

3.1.2 The Scramjet 13 

3.2.3 The Rocket 14 

4. SOLUTION APPROACH 15 

5. RESULTS 16 

5.1 Effect of Acceleration Constraints 16 

5.2 Effect of Available Thrust Level 17 

5.3 Effect of Fuel Specific Impulse 17 

6. CONCLUSIONS AND WORK IN PROGRESS 18 

REFERENCES 19 

TABLES 21 

FIGURES 23 

APPENDIX 33 

A.O INTRODUCTION 33 

A. 1 AMBIENT AND FREE-STREAM CONDITIONS 33 

A. 2 THE AFTERBURNING TURBOJET 33 

A. 2.1 The Intake 34 

A. 2. 2 The Compressor 34 

A. 2. 3 The Combustion Chamber 35 

A. 2.4 The Turbine 35 

A. 2. 5 The Afterburner 36 

A. 2.6 The Nozzle 36 

A. 3 THE RAMJET 37 

A. 4 THE SCRAMJET 38 

A. 4.1 The Intake 38 

A. 4. 2 The Combustion Chamber 38 

A.4.3 The Nozzle 39 

A. 5 CONCLUSION 40 


2 


ACKNOWLEDGEMENTS 

The authors thank the United States Air Force for allowing us to use the OTIS computer 
program. 


NOTATION 


A 

a 

C 

C D 

Cl 

Cy 

S 

ER 

f 

G 

g 

H 

h 

fe 

L 

M 

m 

P 

Px 

Q 

q 

R 

r 

S 

T 

t 

V 

w 

w 

X 


a 

P 

Y 

F 

0 

P 

K 

Subscripts 

a 

c 

KE 

R 

t 

0 


Cross Section Area 
Acceleration 

State and Control Inequality Constraint 

Drag Coefficient 

Lift Coefficient 

Velocity Coefficient 

Constant Pressure Specific Heat 

Drag Force 

Expansion Ratio 

Fuel to Air Ratio 

Acceleration Load Limit 

Gravitational Acceleration 

Heating Value, Hamiltonian 

Altitude 

Fuel Specific Impulse 
Process Efficiency 
Lift Force 
Mach Number 
Mass, Mass Flow 
Pressure 

Adjoint Variable to the State X 
Heating Rate 
Dynamic Pressure 

Earth Equatorial Radius, Gas Constant Air 

Radial Distance to the Center of the Earth, Nose Radius 

Surface Area, State Inequality Constraint 

Thrust, Absolute Temperature 

Time 

Velocity 

Weight 

Jet Velocity 

Down Range 

Angle-of-Attack 

Direction Angle of the Resultant Aerodynamic Force, tan|J = L/D 
Flight-Path-Angle, Ratio of Specific Heats 
Efficiency 
Thrust Pitch Angle 

Gravitational Constant, Lagrange Multiplier 
Linear Throttle 


Ambient 
Combustor 
Kinetic Energy 
Pressure Recovery 
Stagnation Condition 
Standard, Entry 



SUMMARY OF WORK COMPLETED 


The first step in our approach to developing guidance laws for a horizontal take-off, 
air-breathing single-stage-to-orbit vehicle is to characterize the minimum-fuel ascent 
trajectories. We have developed the capability to generate constrained, minimum-fuel ascent 
trajectories for a single-stage-to-orbit vehicle. A key component of this capability is the general 
purpose trajectory optimization program OTIS, developed by the Boeing Company for the Air 
Force. The pre-production version, OTTS 0.96, has been installed and run on a Convex C-l 
at Princeton University. A more recent version, OTIS 1.2, is to be installed in September 
1989. We have developed supporting software for generating initial guesses. By employing a 
direct solution method, state inequality constraints, inherent features of the ascent problem, can 
be routinely handled without resorting to reduced-order approximations for the vehicle point- 
mass dynamics. 

A propulsion model has been developed covering the entire flight envelope of a 
single-stage-to-orbit vehicle. Three separate propulsion modes, corresponding to an 
after-burning turbojet, a ramjet and a scramjet, are used in the air-breathing propulsion phase. 
For each mode, propulsion performance is modeled using a simple one-dimensional flow 
analysis. By varying the propulsion model parameters over a reasonable range, one obtains a 
family of propulsion systems exhibiting certain salient propulsion characteristics covering the 
range of anticipated propulsion system behavior. A variable thrust rocket engine is 
implemented for orbit insertion and flight in the upper regions of the atmosphere. 

The GHAME aerodynamic model of a hypersonic air-breathing single-stage-to-orbit vehicle 
has been obtained and implemented. Note that, though the present aerodynamic and propulsion 
models provide a reasonable starting point for our investigations, other models can be readily 
implemented. 

- Preliminary results pertaining to the effects of variations in acceleration constraints, 
available thrust level and fuel specific impulse on the shape of the minimum-fuel ascent 
trajectories have been obtained. The results show that, if the air-breathing engines are sized for 
acceleration to orbital velocity, it is the acceleration constraint rather than the dynamic pressure 
constraint that is active during ascent. Realistic acceleration constraints can be accommodated 
by proper adjustment of the throttle and angle of attack with only minor mass penalties. The 
acceleration constraint is accommodated by increasing the angle of attack to raise the flight path 
during potentially high acceleration periods and throttling back when the acceleration limit is 
reached. There are finite trajectory segments on the acceleration boundary. The response to 
decreasing the level of available thrust is that the vehicle flies uniformly at lower altitude. The 
flight time is increased, as is the fuel consumption. Changing the I sp has little effect on the 
optimal trajectory and controls, but does increase the fuel consumption as expected. 



1. INTRODUCTION 


Horizontal take-off, single-stage vehicles, using air-breathing as well as rocket propulsion, 
may offer a more economical means of transporting payloads to orbit. The anticipated 
advantages are the operational flexibility of horizontal take-off, the operational simplicity of a 
single-stage, and the propellant mass reduction that results from using air-breathing engines. 
The ultimate objective is to rely completely on the air-breathing engines to accelerate to orbital 
speed. Effective use of air-breathing propulsion over the entire supersonic and hypersonic 
speed range requires a ramjet with both a subsonic-combustion mode and a 
supersonic-combustion mode. To accelerate during the subsonic and transonic phases some 
additional propulsion mode(s) is necessary. Further discussion of the propulsion system 
development can be found in [1]. The effective use of air-breathing engines leads to lower 
altitude ascent trajectories as compared to those of the Space Shuttles. Higher dynamic 
pressure, higher surface temperatures, and higher acceleration loads result. 

A number of studies have been conducted to determine the performance potential of an 
air-breathing SSTO vehicle. The study by Schoettle [2] provides useful background for our 
work. Schoettle compared rocket versus combined rocket and ramjet/scramjet performance and 
concluded that the use of a ramjet engine(s) allows substantial propellant savings, but is heavily 
penalized by the added weight of the ramjet engine and the structural reinforcements required to 
accommodate the higher dynamic pressures. The maximum payload delivery trajectory 
involved dynamic pressures of up to 5850 psf. Only by accommodating dynamic pressures up 
to at least 1044 psf could a payload delivery advantage over a rocket driven vehicle be realized. 
As a point of reference, the Space Shuttles do not experience dynamic pressures above 650 psf 
during their ascent. The payload mass versus peak dynamic pressure curve had a steep slope 
up to 1566 psf, at which point the payload mass was 80% of the maximum achievable value. 
Between 1566 psf and 5850 psf, the slope was small and relatively constant. The results of 
this study by Schoettle and other similar studies offer incentive for continuing research towards 
the development of an air-breathing SSTO vehicle. 

Our objective is to develop guidance logic for a horizontal take-off, air-breathing SSTO 
vehicle. As a first step, we aim to characterize the minimum-fuel trajectories that an air- 
breathing SSTO vehicle should ideally fly to achieve maximum performance. Because the 
optimal trajectory characteristics are critically dependent on the propulsion system and 
aerodynamic models one adopts, and because there is significant uncertainty as to what these 
models should be, we will parameterize our ignorance as much as possible, so as to explore the 
range of possibilities and determine the sensitivities to various features of the models. We will 
also explore the effects of dynamic pressure, heating, and acceleration constraints on the 
character of the optimal trajectories. Our guidance approach will then follow naturally from an 
understanding of the ideal optimal trajectories and controls. That is, the characteristics and 
appropriate approximations on which to base the guidance scheme will be inferred from the 
results. 

In this report, we lay the groundwork for our characterization of minimum-fuel 
trajectories for an air-breathing SSTO vehicle. The problem formulation is given. The 
propulsion system and aerodynamic models are discussed. The solution approach is outlined. 
And some preliminary results are presented. 


2. PROBLEM FORMULATION 


In all its generality, the determination of an ascent trajectory, let alone an optimal one, for 
any launch vehicle constitutes an elaborate problem in itself. With the added complexity of 
sustained flight in the atmosphere for an extended period of time, inherent to an air-breathing 
SSTO vehicle, the problem becomes even more formidable. The purpose of this research 
however, is to obtain a characterization of minimum-fuel ascent trajectories. The optimal 
trajectory characteristics will provide the basis for the subsequent development of a guidance 
scheme. Therefore, the introduction of some simplifying assumptions which would reduce the 
overall complexity of the ascent problem without significantly affecting its character is justified. 

First, the Earth is assumed to be modeled as a homogeneous sphere. Neglecting the 
oblateness and inhomogeneity of the Earth implies a central gravity field and allows us to use 
the simple inverse square law (g = |i/r 2 ) for the gravitational acceleration. 

Second, the Earth is assumed to be non-rotating. This simplification is acceptable 
considering the very small accelerations experienced due to the Coriolis force. Furthermore, if 
this force was included the orientation of the ascent trajectory with respect to the Earth 
rotational axis would come into play, complicating a general trajectory characterization. 

Third, the atmosphere is taken to be stationary and in what follows the use of a standard 
atmosphere providing ambient conditions as a function of altitude only has been assumed. 

2.1 Equations of Motion 

At this stage we are not considering ascent from a specific launch site to a specific target 
orbit, but the simpler problem of ascending to a generic low Earth orbit. Intuitively, the 
minimum-fuel ascent trajectory for this case will not exhibit out-of-plane motion, that is, the 
ascent trajectory will be confined to a great circle plane and the heading is fixed. Imposing this 
restriction a priori eliminates the bank angle as a control variable, simplifying the control 
problem even further. Subsequent work dealing with specific launch site-target orbit pairings 
will need to take out-of-plane motion into consideration. With this additional simplification and 
the usual no-slip condition, the dynamic equations describing the motion of the vehicle (or 
more precise the motion of the center of mass of the vehicle) are [3,4] 
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where R is the radius of the Earth and go = lf/R 2 . Observe that the direction of the thnist is 
fixed with respect to the body axes and lies in the vehicle’s plane of symmetry. The line of 


thrust is offset with a fixed thrust pitch angle 0 with respect to the line of zero angle-of-attack. 
The aerodynamic forces are described by 

L = Cl q S re f D = Cd q S re f (6) 

where the lift coefficient (Cl) and the drag coefficient (Cd) are functions of the angle-of-attack 
(a) and Mach number (M), q is the dynamic pressure and S re f is a reference area. The 
propulsive force is given most generally by 


T — ^ j^T max (M, P a » ^ a’ °0 

k k (7a) 

or, equivalently in a standard atmosphere (e.g. the U.S. Standard Atmosphere, 1976 [5]) 
where the ambient pressure and temperature are functions of the altitude (h = r - R) only, by 


T — ^ 7t k T max (M, h, a) 
k k 


(7b) 


where 7Ck is a continuous linear throttle factor ranging from zero to one. The thrust dependence 
on a reflects the influence of incidence on the engine intake performance. The index k allows 
for more than one particular propulsion mode to be effective at the same time: for example, the 
use of a throttleable rocket along with a scramjet in the higher regions of the atmosphere. 

Likewise, the specific impulse (I sp ) is a function of Mach number and either the ambient 
conditions or the altitude in a standard atmosphere 
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The common assumption is made here that the value of the I sp is not affected by throttling of 
the engine. This assumption is not always justified and the I sp can depend in a nonlinear 
fashion on the throttle setting. However, determining I S p as a function of the linear throttle 7t 
involves the non-trivial inversion of the propulsion model and seems unwarranted at this stage. 

The vehicle dynamics are thus controlled by the angle-of-attack a and the set of linear 
throttles Ttk- 


2.2 Initial and Final Conditions 


Ideally one would like to consider the ascent from release of the brakes at the beginning of 
the runway, to insertion into low Earth orbit (LEO). In practice however the execution of a 
minimum fuel ascent will be restricted by safety and operational considerations to commence 
after take-off and initial-climb out. Additionally, the acceleration to supersonic speeds 
employing a turbojet-like mode of propulsion is well understood and this phase of the ascent to 
orbit accounts for only a small portion of the total fuel consumption. Hence our work has been 
focused on the high-speed phase of the ascent. 

Similarly, the final insertion into low Earth orbit using a rocket orbital maneuvering system 
is well understood and actual determination of the orbit insertion might be computationally 
intensive due to extreme sensitivity of various orbital parameters to the insertion parameters. 
Instead, acceleration to circular speed at the edge of the sensible atmosphere is used to 
benchmark the various optimal ascent trajectories. 



Subsequent work may include a more detailed look at the coupling between the trunk ascent 
trajectory and the initial climb-out and orbital insertion phase respectively. 

2.3 Performance Index 

Since this work is concerned with minimum-fuel ascent to orbit the obvious performance 
index to be minimized is the fuel expenditure. For a vehicle with a fixed gross take-off weight 
this is equivalent to maximizing final weight. The ultimate purpose of the air-breathing ascent is 
to maximize payload/fuel mass fraction. To this end one should make the gross take-off weight 
a design (i.e. free) variable. This has not been taken up here since, if done properly, it would 
amount to a full-fledged design study. Our study is limited to the case of a fixed gross take-off 
weight and the minimum-fuel ascent trajectories are obtained by formally maximizing the final 
mass. 

2.4 State and Control Constraints 

Vehicle and crew considerations dictate certain constraints that must be placed on the states 
and/or controls. The throttle settings must take on values between 0 and 1. The angle of attack 
must take on values between a m i n and a max . The value of a max corresponds to either the stall 
boundary or a propulsion system dictated boundary. The minimum angle of attack is dictated 
by the propulsion system characteristics or the limits of our aerodynamic data. <x m in is typically 
some small negative angle. 

We consider two state and control inequality constraints that place limits on the axial and 
normal acceleration loads 

Ci = | 1/m (rtTmax — (L 2 + D 2 )!/ 2 cos(oc + [})) | - Gaxial — 0 (9) 

C2= | 1/m ((L 2 + D 2 ) 1 / 2 sin(a + P)) | - G n0 nnal ^ 0 (10) 

where tan (3 = L/D and G ax ial and G n0 rmal are the imposed acceleration limits. A less specific 
option for constraining the acceleration load is to limit the total acceleration load 

C = 1/m (tt 2 T max 2 + L 2 + D 2 + 2rcT max (Lsina - Dcosa)) 1 / 2 - Gtotal ^ 0 

We also consider the two first-order state inequality constraints that place limits on the 
heating rate and dynamic pressure 

Si= Kpl/ 2 V3-Q max £ 0 (11) 

S 2 = \ p 1 / 2 V 2 - q m ax < 0 (12) 

Clearly, the set of constraints could be greatly expanded, but this would require specific 
vehicle knowledge. In order to keep our study generic and not over-complicate the 
computations, we will not go beyond the above set of constraints. 

2.5 Optimality Conditions 

Even though we will be employing a direct method of solution, some insight into the 
characteristics of the minimum-fuel trajectories and controls may be obtained from deriving the 
optimality conditions that would form the basis of an indirect method of solution. 



The Hamiltonian for the constrained minimum-fuel ascent problem is 
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The equations for the adjoint variables are 
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If the optimal trajectory contains an arc that is on the boundary of one of the first-order state 
constraints, conditions of the form [6] 
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apply at the entry or exit from the boundary arc, where tk is the entry or exit time and X r and 
X v are scalar Lagrange multipliers. The entry and exit points may or may not be comers, i.e., 
places where some or all of the controls are discontinuous. 

For flight away from the constraint boundaries, the optimal controls are as follows. Based 
on the Maximum Principle, the values of the controls at each point along the optimal trajectory 



should maximize H. Taking the partial derivatives of H with respect to the controls and 
expressing the result in terms of the state variables used in Eqs. (l)-(3) and (5) and the 
corresponding adjoint variables, we obtain 
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Since the control n appears linearly in H and 3H/37I is its coefficient, H is maximized by 
choosing K according to the rule 


k = 0 if 3H/a7t<0 

7t = 1 if 3H/a7t>0 


(24) 


If 3H/57t = 0 over a finite interval of time, the control is singular and must be determined 
from higher order conditions. We will show below that the singular case can be ruled out. The 
control a is bounded from above and below. H may be maximized by <x m in> «max> or an 
interior value of a that satisfies the implicit equation 
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The Hessian of H with respect to the controls, denoted by H uu , has components 
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If an interior value of a satisfying Eq. (25) maximizes H, the second partial given in 
Eq. (28) should be less than or equal to zero; if it is equal to zero, higher derivatives with 
respect to a would have to be checked. 

In order to determine whether a singular thrust arc can be optimal, we examine the 
determinant of H uu . From Eqs. (26)-(28), it follows that 
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Thus det(H uu ) < 0. If det(H uu ) < 0, H uu is indefinite and the singular arc is not optimal. If 
det(H uu ) = 0, then the singular arc may be optimal and additional conditions would have to be 
checked. The first factor on the right-hand-side of Eq. (29) is never zero based on physical 
reasoning; the second factor is zero only if p v = Py = 0 or 


tana - (30) 

Vp v 

The first condition implies that p r = p m = 0 as well and is inconsistent with an optimal 
solution. Referring to Eq. (25), the second condition can only be satisfied in the limit as the 
density p (and hence q) goes to zero, i.e., outside the atmosphere. Consequently, we conclude 
that singular thrust arcs are not optimal during the ascent through the atmosphere. Corban, 
Calise, and Flandro[28] reached the same solution for a reduced-order model of the point mass 
dynamics. 

Transversality conditions, integrals of the motion, and the optimal controls for flight on the 
constraint boundaries will not be developed here, since no further insight is expected. We only 
remark that intermediate values of throttle are possible along constraint boundaries. The results 
to be described later show that indeed intermediate values of throttle are used to fly along the 
boundary of the acceleration constraint. 

In the above development, the issue of choosing the optimal propulsion mode has not been 
addressed. Accounting for the existence of different propulsion modes, we rewrite the 
Hamiltonian as 
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where k is the number of propulsion modes and the subscript "reduced" denotes that the terms 
not involved in determining the optimal throttle setting have been suppressed. In the second 
line of the equation, the reduced Hamiltonian is expressed in the form of an inner product 
(actually there are k inner products). At this point in our investigations, we are assuming that 
only one propulsion mode is active at any given time. The optimal mode is the one associated 
with the largest of the k inner products. Since the angle of attack appears in these inner 
products and the k throttle settings appear in the equations for determining the optimal angle of 
attack, the optimal controls must be determined by simultaneous maximization of the 
Hamiltonian. At points along the optimal trajectory where there is a switch from one propulsion 
mode to another, there will in general be comers where there are discontinuities in the time 
derivatives of the speed and flight path angle. 


3. VEHICLE MODELING 


The current vehicle configuration is conform a series of NASA reports on a Hypersonic 
Research Airplane Concept [7-10]. The airframe is expected to be a cross between a flying 
body and a delta wing with a belly-mounted multi-mode engine. The gross take-off weight of 
the SSTO vehicle is estimated to be 300,000 lbf of which approximately 60 % is fuel [11]. The 
overall length of the vehicle is approximately 233.4 ft and its aerodynamic reference area is 

6.000 sqft. 

Liquid hydrogen is used as fuel over the entire speed range since it allows for the high 
effective specific impulse that makes air-breathing propulsion a viable option and is suitable for 
supersonic combustion. Recently, slush hydrogen has been mentioned as an operationally 
more attractive form of the hydrogen fuel [12]. 

3.1 The Aerodynamic Characteristics 

Currently, there is a dearth of up-to-date aerodynamic data concerning lifting vehicles 
flying at hypersonic speeds. However, representative values for the lift and drag coefficients 
covering the entire flight regime of a SSTO vehicle, have been obtained from . The generic 
hypersonic aerodynamic data therefore contains certain realistic aerodynamic anomalies. 

The present study is concerned with guidance for performance, that is, it considers the 
motion of the vehicle's center of mass. The Cl and Cd values used in this point mass model 
are the trimmed (equilibrium) values as functions of Mach number and angle-of-attack a; the 
rotational dynamics are not accounted for.ones obtained by instantaneous trimming of the 
vehicle. We currently employ the trim map found in the Generic Hypersonic Aerodynamic 
Model Example (GHAME) [11] which covers the entire flight regime of a SSTO vehicle. This 
model is an amalgamation of the characteristics of the space shuttles, the X-24C lifting body 
and a 70° swept double delta wing configuration. In Figure 1, Cl and Cd are plotted as a 
functions of Mach number for a selected number of angles-of-attack. These plots are obtained 
from the original data grid by quadratic interpolation over a for fixed Mach number and 
subsequent determination of cubic splines for the quadratic coefficients over the Mach number 
range. Actual numerical calculations however employ local 2-dimensional quintic interpolation 
applied to the original data grid. 

3.2 The Propulsion System 

The aero-engine for an air-breathing SSTO vehicle comprises several distinct propulsion 
modes, each suited to operate over a specific speed range. Our present propulsion model 
consists of an afterburning single-flow turbojet, a ramjet and a scramjet. The model is 
designed to provide maximum specific thrust and the corresponding fuel specific impulse. The 
respective modes are analyzed using simple one-dimensional flow analysis for perfect gases, 
employing appropriately constant efficiencies [13-18]. 

A simple turbojet model has been developed for completeness and may or may not 
accurately represent the actual low speed propulsion mode; there is still much uncertainty about 
the kind of propulsive device to be used for take-off and climb-out. 

The emphasis of our work has been on the high-speed phase of the ascent after initial 
acceleration to supersonic speeds. At these speeds a ramjet with either subsonic or supersonic 
combustion is most suitable. Though we do not want to concern ourselves with the mechanical 
particulars of such an engine, one may visualize an integrated device where upstream transfer 
of the combustion process effects a switch from a ramjet with conventional subsonic 
combustion to a ramjet with supersonic combustion; a scramjet. A brief description of the 
respective ramjet propulsion modes is presented in what follows, for more detail the reader is 
referred to the Appendix. 
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The advantage of the present propulsion model based on fundamental thermodynamic 
relations is that it allows us to obtain different propulsion characteristics by changing specific 
engine efficiencies and parameters. 

At high speeds the entire vehicle becomes a flying engine, its forebody functioning as a 
compression surface and its trailing edge as an expansion surface. The inlet and exit 
characteristics of the propulsion device can therefore not be modeled accurately by a simplified 
analysis or existing design methods (e.g. additive drag, intake efficiency and nozzle losses). 
However, it is reasonable to expect that the variable geometry inlet is designed to reduce the 
effect of inlet spillage on the performance and that the intake is always matched to the engine. 
Hence, for the moment, additive drag has been neglected. Similarly, the 3D-nozzle is 
represented only by its thermodynamic workings. A simple estimate for the internal friction of 
the engine has been incorporated. 

The effect of incidence on the engine performance with belly-mounted engine pods is 
twofold [19]. First, increased incidence leads to a larger mass capture by forebody-intake 
system and thus will enhance the thrust. On the other hand increased incidence produces a 
stronger bow shock and will decrease the efficiency of the intake process thus decreasing the 
thrust. At very high incidence the loss in efficiency will outweigh the gain in massflow since 
the latter has an upperbound dictated by the engine size. The two opposite trends will therefore 
most likely result in either an immediate roll-off in thrust or an initial increase followed by a 
roll-off. 

3.1.1 The Ramjet. The airflow is diffused adiabatically to subsonic speeds with constant 
kinetic efficiency. The airflow enters the combustor at constant Mach number. Combustion 
takes place with constant burning efficiency and a fixed drop in stagnation pressure. The 
amount of fuel is determined by the maximum combustor exit temperature. 

Next the airflow is accelerated again to supersonic speeds through a continuously 
adjustable convergent-divergent nozzle. Ideally the flow should expand to ambient pressure for 
maximum performance, however, the nozzle's expansion ratio is limited by the maximum size 
of the exit area. 

The specific thrust and specific impulse determined from this model for a couple of 
altitudes in the 1976 U.S. Standard Atmosphere are presented in Figure 2, for the specific set 
of ramjet parameter values presented in Table 1. Towards Mach = 7 the combustor stagnation 
entry temperature approaches the fixed combustor stagnation exit temperature, resulting in a 
rapid decline in thrust and a steep drop in specific impulse. Observe furthermore the relative 
insensitivity of the specific impulse as a function of Mach number to changes in altitude as 
opposed to the thrust. 

3.1.2 The Scramjet. In this propulsion mode the airflow is diffused to supersonic speeds 
with constant process efficiency. Swithenbank [14] recommends a diffusion ratio of 3.0 for 
optimal specific impulse, hence the intake delivery Mach number is taken to be one third of the 
flight Mach number. 

In our original scramjet model, the supersonic combustion takes place at constant 
combustor cross section area. An approximate optimal fuel-to-air ratio as a function of flight 
Mach number is obtained from [14]. The fuel mixture is weaker than stochiometric at flight 
speeds below M = 10 and richer at speeds approaching circular velocity. 

m f / m . = 0. 003 xM M<12. 5 

fuel air 



M >12.5 


( 32 ) 


In case this fuel-to-air ratio leads to thermal choking the amount of fuel is reduced to the critical 
amount, for which the flow reaches marginal choking conditions at the combustor exit. 
Thermal choking occurs for flight Mach numbers less than approximately 10.7. This resulted 
in an inadequate level of thrust in the lower Mach number range of M = 4 to 6. The range 
where, according to reports [20], the scramjet is expected to kick in. In our revised model, the 
supersonic combustion can take place in a variable cross section area burner. The additional 
parameter also allows for the limiting case of constant pressure combustion, which is thought 
to be more conducive to the combustion process. 

Since the flow is supersonic, the nozzle consists of a divergent section only. The expansion 
ratio is again fixed by the size of the exit area. The nozzle efficiency as expressed by the value 
of the velocity coefficient C v turns out to have a significant effect on the form of the propulsion 
characteristic. Over a relatively narrow range of C v values the specific thrust characteristic is 
found to change from monotonically increasing with Mach number to rolling off at high Mach 
numbers. Since there is still much uncertainty about the true behavior of the scramjet at very 
high speeds, it is important to determine the effect of this variation in thrust character on the 
make-up of the ascent trajectory. 

The specific thrust and specific impulse determined from this model for a couple of 
altitudes in the U.S. Standard Atmosphere, 1976 are presented in Figure 3 for the original 
model and the specific values for the scramjet parameters presented in Table 2. Figure 4 
shows the characteristics of the revised model at two different values of the nozzle velocity 
coefficient C v , with the specific parameter values presented in Table 3. 

The airflow through the combustor is thermally choked below Mach = 10.7 in the original 
model, this shows up as a corner in the thrust and I sp graphs. With the revised model this 
comer indicates the transition from a variable cross-sectional area combustor to a constant 
cross-sectional area one. The second comer in both graphs is related to the slope discontinuity 
in the fuel-to-air ratio at M = 12.5. The rise in the specific impulse of the scramjet at high 
Mach numbers is caused by the specific impulse contribution due to the kinetic energy of the 
fuel. This becomes a significant contribution at high speeds. Furthermore, the revised model is 
seen to predict I sp levels of 3200 to 4000 seconds, consistent with those cited in [21]. 

3.2.3 The Rocket. In order to be able to maneuver in space the SSTO vehicle must have a 
rocket engine. Depending on the performance of the scramjet at high speeds and high altitudes 
there might also be a need to employ a rocket before exiting the atmosphere. 

The RL-10 rocket engine for the Centaur upper stage is suggested by Calise et al [22] to be 
representative of a rocket engine to be installed in a SSTO vehicle. It is rated at 15,000 lbf 
thrust and a I sp of 444 sec at 200,000 ft [23]. Back pressure will be taken into account 
employing an exhaust area of 11.5 sqft. 


4. SOLUTION APPROACH 


Flight path optimization, based on a calculus of variations formulation, entails the solution 
of a nonlinear two point boundary value problem [6]. This is a difficult task which is 
exacerbated when state and control inequality constraints are imposed. Because of the critical 
role of state and control inequality constraints in the present flight path optimization problem, 
we sought an optimization algorithm that could readily accommodate such constraints. We are 
currendy using a computer program entitled Optimal Trajectories by Implicit Simulation (OTIS) 
[24], installed on a Convex C-l. Thus far it has served us well. 

The approach, on which OTIS is based, is to formulate the optimal control problem as a 
nonlinear programming problem using implicit integration of the trajectory by a collocation 
method. The states along the trajectory as well as the controls are represented by cubic splines. 
The collection of breakpoints (or nodes) of the cubic segments constitutes a discretization of the 
trajectory. By adjoining the defects of the collocation at the midpoint between successive nodes 
to the original performance index, one obtains an augmented performance index, subject to the 
original boundary conditions and constraints evaluated at the nodes and the midpoints. In this 
manner the original optimal control problem is stated in the form of a constrained nonlinear 
programming problem [25]. The independent variables are the states and the controls at the 
nodes plus the length of the time interval. To solve such a programming problem a host of 
methods is available, see [26]. The present version of OTIS employs the NPSOL 2.1 package 
[27]. 

OTIS needs an initial estimate, preferably a feasible trajectory. We obtained our original 
feasible trajectory by piecing together segments of flight at constant dynamic pressure, constant 
altitude and constant throttle setting. The trajectories of previous runs were used for subsequent 
initial estimates. 

Observe that the augmented performance index is only locally minimized and depends on 
the initial estimate, the scaling of variables and constraints and the specific node distribution. In 
order to accept any result as a locally optimal solution, one needs to verify relative insensitivity 
to the node distribution. The total number of nodes employed currently is 65, which is pushing 
the capacity of our machine. 


5. RESULTS 


We present some preliminary results which only begin our characterization of the 
minimum-fuel ascent trajectories. The conditions at the end of the initial climb-out and 
acceleration phase, the phase in which the turbojet propulsion mode is assumed to be used, are 
estimated to be Mach 2 in speed, 1,000 psf in dynamic pressure, and zero degrees in flight 
path angle. The weight of the vehicle is estimated to be reduced to 290,000 lbf at the end of the 
climb-out, from 300,000 lbf gross take-off weight. These quantities determine the initial 
conditions for the minimum-fuel ascents (see Table 4a). The final conditions correspond to a 
circular orbit at the edge of the sensible atmosphere, which is thought to occur at 259,000 ft, 
leading to a final speed and flight-path-angle of 25,777 fps and 0° respectively (see Table 4b). 

In order not to limit, a priori the vehicle’s flight envelope, the engines are sized so that 
orbital insertion using air-breathing propulsion only is feasible. With insertion into a 50 nm by 
100 nm transfer orbit as a typical target orbit, a combustor cross-sectional area of 400 ft 2 was 
found to provide an adequate level of thrust of 4,000 lbf at the insertion point with the original 
propulsion model. Together with data mentioned in earlier chapters, this defines the 
(preliminary) baseline vehicle (see Table4c). Additionally, the thrust level does not depend on 
the angle-of-attack, i.e., the sensitivity of thrust with respect to incidence effects has not been 
considered. The thrust pitch angle is identically zero. 

In the numerical computations the angle-of-attack range was restricted to the range of our 
aerodynamic data (-3° to 21°) and the a rate of change was limited as to require one second to 
traverse the entire a range. Likewise, the throttle rate of change was restricted to be less than 
20 %/sec, precluding the occurrence of bang-bang and chattering controls. 

5.1 Effect of Acceleration Constraints 

The baseline model with engines sized for orbital insertion and capable of generating high 
levels of thrust in the lower regions of the atmosphere will consequently be able to accelerate 
very rapidly. However, as mentioned in section 2.4, there are limits to the allowable level of 
sustained acceleration. It is therefore of interest to assess the effect of the magnitude of the 
acceleration limit on the character of the ascent trajectory and to evaluate the relative penalty 
incurred using a more stringent limit. Three acceleration limits where implemented varying 
from 5g’s to 3g’s. 

The results in Table 5a show that, as expected, the performance as measured by the value 
of the final mass diminishes as the acceleration constraint becomes more stringent, while the 
time of ascent increases. The penalty one has to pay for the convenience of ascending to orbit at 
3g's as opposed to 5g's is however a meager 1,000 lbf. 

Additionally, it is interesting that the Mach number and altitude at which the switch is made 
from ramjet to scramjet shows little sensitivity to the value of a m ax- This is also evident from 
Figure 5. Referring back to Figures 1 and 2, note that the switch occurs in the Mach number 
range where the I sp of the ramjet starts falling off but is still superior to the scramjet I™ which 
is near its maximum value. The bottom curve in Figure 5 represents a constant dynamic 
pressure of 1,500 psf. It shows that along all three ascent trajectories the dynamic pressure is 
less than 1,500 psf, although no dynamic pressure constraint has been imposed. 

Figures 6a and b show another general characteristic: the vehicle uses full throttle unless it 
has reached the acceleration constraint, in which case it throttles back in order to follow the 
constraint. The same figures also illustrate the switch between propulsion modes. At both sides 
of the switch the propulsion system is at full throttle. The switch occurs apparently when the 
drop in I sp in going from ramjet to scramjet is outweighed by the higher scramjet thrust. 

Figure 7 presents the other control variable, the angle-of-attack. The character of the 
angle-of-attack vs. Mach number curve is similar for the three acceleration limits. The 


angle-of-attack is generally higher, the lower the acceleration limit. Higher angle-of-attack 
produces higher lift and consequently the higher altitude flight path seen in Figure 5. 

5.2 Effect of Available Thrust Level 

The combustor cross-sectional area for the baseline model roughly corresponds to a frontal 
intake area of 670 sqft [15]. The intake to reference area ratio of 0. 1 1 2 contrasts sharply with 
the area ratio of 0.0195 for the Hypersonic Research Airplane Concept of [7-10]. Although it 
should be observed that the latter is a hypersonic cruise vehicle and not a SSTO boost vehicle, 
the discrepancy points out that the actual engine size could be determined by factors other than 
air-breathing orbital insertion capability. Therefore, the available thrust level may not be quite 
as high as assumed for the baseline model. Using the 300 sqft frontal intake area suggested in 
[1 1], a combustor cross-sectional area of 180 sqft is obtained. The effect of the overall thrust 
level on the minimum-fuel ascent is estimated, evaluating three engine sizes varying from the 
baseline 400 sqft combustor to a 180 sqft one. A 3g’s acceleration limit is employed in all 
three cases. 

The results in Table 5b show a 33% increase in final time, resulting in increased fuel 
consumption. The propulsion mode switching still occurs in the same Mach number range, 
while the linear throttle exhibits the same kind of generic behavior as noted in the previous 
section, see Figure 6c. However, as illustrated by Figure 8, the vehicle with the reduced 
engine size shows an overall tendency to fly lower at higher levels of dynamic pressure. This 
tendency most likely stems from the fact that thrust is proportional to dynamic pressure. 

5.3 Effect of Fuel Specific Impulse 

Presently, there is much uncertainty concerning the performance of the sub/supersonic 
combustion ramjet. Therefore, it is of interest to determine what effect a drop in overall engine 
performance, as expressed by the fuel specific impulse, has on the character of the 
minimum-fuel ascent. Employing the baseline model and again a 3g acceleration limit, three 
trajectories have been determined with I sp levels ranging from 100% to 60% of our baseline 
values. 

The results presented in Table 5c indicate a modest sensitivity to the engine performance as 
far as trajectory make-up is concerned, with a slight increase in altitude at lower I sp (see 
Figure 9), where lower ambient temperatures lead to improved performance. In contrast, the 
effect on the final mass is very profound as was to be expected. The linear throttle still shows 
the same generic behavior, effectively functioning as a direct acceleration control. 


6. CONCLUSIONS AND WORK IN PROGRESS 


We have developed the capability to generate constrained, minimum-fuel ascent trajectories 
for a single- stage-to-orbit vehicle. By employing a direct solution method, state inequality 
constraints, inherent features of the ascent problem, can be routinely handled without resorting 
to reduced-order approximations for the vehicle point-mass dynamics. A model, based on one- 
dimensional flow analysis, has been developed for a multi-mode propulsion system capable of 
turbojet, ramjet, scramjet, and rocket performance. The free parameters of the propulsion 
model can be adjusted to obtain a range of performance. Aerodynamic parameters have been 
obtained from the GHAME. The developed capability enables us to characterize the minimum- 
fuel ascent trajectories. Subsequent guidance law development will be based on the results of 
the characterization. 

Results pertaining to the effects of variations in acceleration constraints, available thrust 
level and fuel specific impulse on the shape of the minimum-fuel ascent trajectories have been 
obtained. The results show that, if the airbreathing engines are sized for acceleration to orbital 
velocity, it is the acceleration constraint rather than the dynamic pressure constraint that is 
active during ascent. Realistic acceleration constraints can be accommodated by proper 
adjustment of the throttle and angle of attack with only minor mass penalties. The acceleration 
constraint is accommodated by increasing the angle of attack to raise the flight path during 
potentially high acceleration periods and throttling back when the acceleration limit is reached. 
There are finite trajectory segments on the acceleration boundary. The response to decreasing 
the level of available thrust is that the vehicle flies uniformly at lower altitude; the flight time is 
increased, as is the fuel consumption. Changing the I S p has little effect on the optimal trajectory 
and controls, but does increase the fuel consumption as expected. 

Our initial minimum-fuel trajectory results showing acceleration to be the limiting constraint 
are inconsistent with those of other investigators [2,22,28] who have found dynamic pressure 
to be the limiting constraint. Since the work reported in Chapter 5 was completed, a revised 
scramjet model has been developed. This revised model allows us to investigate the effect of 
thrust roll-off at high speeds on the characteristics of the minimum-fuel ascent. We are 
currently computing minimum-fuel trajectories, for which the rocket engine is allowed to 
augment the air-breathing thrust capability during the ascent. With this approach, we avoid 
having to size the air-breathing engines for acceleration to orbital velocity. As discussed in 
Section 5.2, lower air-breathing engine thrust levels lead to higher dynamic pressures. Thus 
we should obtain results that are more in line with those of other investigators. 

We also intend to investigate the effects of making the thrust level dependent on 
angle-of-attack or the thrust pitch angle. Scaling the thrust by a quadratic function of angle of 
attack will allow us to investigate several types of thrust variation by using different values of 
the coefficients in the quadratic function. For the proposed type of air-breathing SSTO vehicles 
where the fore- and aftbody of the vehicle constitute integral parts of the propulsion system, the 
direction of thrust will vary with flight conditions. We plan to develop a means simulating this 
feature and then to determine the consequences. 

Much useful insight into the characteristics of the minimum-fuel ascent trajectories has been 
gained using reduced-order models [22,28]. An important question that needs to be answered 
is: How well do the optimal trajectories for the reduced-order models approximate the optimal 
trajectories for the full-order model? We intend to answer this question and are proceeding to 
develop software to generate reduced-order solutions. 

Finally, we plan to spend much of the next year on initial guidance law development. 
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Table 1 Ramjet Mode Engine Parameters 


Intake Kinematic Efficiency 

0.94 

Combustor Entry Mach Number 

0.5 

Combustion Efficiency 

0.96 

Combustor Stagnation Pressure Loss 

3 % 

Combustor Stagnation Exit Temperature 

2500 K 

Nozzle Velocity Coefficient 

0.97 

Nozzle Exit- Combustor Cross Section Area Ratio 

6 

Table 2 Scramjet Mode Engine Parameters (Original Model) 

Intake Process Efficiency 

0.93 

Diffusion Ratio M0/M3 

3 

Combustion Efficiency 

0.96 

Nozzle Velocity Coefficient 

0.96 

Nozzle Exit-Combustor Cross Section Area Ratio 

6 

Table 3 Scramjet Mode Engine Parameters (Revised Model) 

Intake Process Efficiency 

0.92 

Diffusion Ratio M0/M3 

3 

Combustor Area Change Parameter 

E= 0.0 M < 3.0 

E =(M-3)/7.7 3.0 < M < 10.7 
E = 1.0 M > 10.7 

Combustion Efficiency 

0.93 

Nozzle Velocity Coefficient 

0.93/0.91 

Nozzle Exit-Combustor Cross Section Area Ratio 

6 

Table 4a Initial Conditions 

Mass 

9013.5 slugs 

Mach Number 

2.0 

Dynamic Pressure 

1,000 lbf/ft 2 

Table 4b Final Conditions 

Altitude 

259,000 ft 

Velocity 

25,777 ft/sec 

Flight Path Angle 

0.0° 


Table 4c Baseline Vehicle Model 


Gross Take-Off Weight 

300,000 lbf 

Reference Area 

6,000 ft 2 

Combustor Cross Section Area 

400 ft 2 

Fuel Specific Impulse 

100% 
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Table 5a Effect of Acceleration Constraint (Baseline Model) 


amax 

Mswitch 

hswitch 

9peak 

mfinal 

tfinal 

(g) 


(ft) 

(lbf/ft 2 ) 

(slugs) 

(sec) 

3.0 

6.59 

100,987 

1,025 

5,844 

557.1 

4.0 

6.63 

98,902 

1,134 

5,863 

524.7 

5.0 

6.66 

97,373 

1,273 

5,873 

504.5 

Table 5b Effect of Available Thrust Level (amax = 3 g's) 



Acomb 

Mswitch 

hswitch 

Qpeak 

mfinal 

tfinal 

(ft 2 ) 


(ft) 

(lbf/ft 2 ) 

(slugs) 

(sec) 

400 

6.59 

100,987 

1,025 

5,844 

557.1 

275 

6.60 

94,936 

1,087 

5,776 

614.6 

180 

6.59 

84,252 

1,650 

5,641 

739.4 

Table 5c Effect of Specific Impulse (amax = 

3 g's) 



Isp 

^switch 

hswitch 

9peak 

mfinal 

tfinal 

(%) 


(ft) 

(lbf/ft 2 ) 

(slugs) 

(sec) 

100 

6.59 

100,987 

1,025 

5,844 

557.1 

80 

6.59 

101,045 

1,028 

5,245 

549.7 

60 

6.60 

101,240 

1,029 

4,382 

539.1 
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Figure 9 Effect of Fuel Specific Impulse on the Altitude Profile 
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APPENDIX 


A.O INTRODUCTION 

The purpose of this propulsion model is to render the maximum thrust (T) and specific 
impulse (I S p) as a function of Mach number (M) and altitude (h) in a standard atmosphere, 
covering the entire flight envelope of an air-breathing SSTO vehicle. 

The aero-engine for such a vehicle comprises several distinct propulsion modes, each 
suited to operate over a specified speed range. The propulsion model consists of an 
afterburning single-flow turbojet, a ramjet and a scramjet. The mechanical details of how these 
modes are to be integrated into one device will not be discussed in what follows. It is tacitly 
assumed that such an engine can be construed. 

The respective modes are analyzed using simple one-dimensional flow analysis for perfect 
gasses, employing constant efficiencies throughout the engine. The key however, is to use the 
most appropriate ones, i.e. to use those efficiencies that remain more or less constant over the 
operating range. 

A.l AMBIENT AND FREE-STREAM CONDITIONS 


The ambient conditions, i.e. pressure (p a ), temperature (T a ) and speed of sound (a), are 
completely determined by the altitude in a standard atmosphere. In the current implementation, 
the U.S. Standard Atmosphere, 1976 is being used [4], 

The ambient conditions together with the free-stream Mach number determine the 
free-stream stagnation pressure (pto) and stagnation temperature (Tto) 


p t0 = p 0 V 
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A.2 THE AFTERBURNING TURBOJET 

The basic components of a afterburning turbojet (Figure Al) consist of an intake, 
compressor, combustion chamber or burner, turbine, afterburner and a nozzle. 

In the performance assessment of a turbojet one needs to make a distinction between the 
operation of the engine at its design point and away from it. Since the SSTO vehicle is 
accelerating to Earth-orbit, it operates over a wide range of velocities and altitudes. The off- 
design operation of the engine is therefore its "true" operating condition and should be 
satisfactorily modeled. 

First, the design point operation is analyzed for fixed design parameters such as 
compressor pressure ratio (e c des)> turbine entry temperature (TET) and nozzle expansion ratio. 
Subsequently, using some realistic assumptions about the turbojet operation, estimates can be 
found for the off-design compressor pressure ratio (e c ) and ditto mass flow (m), which are 
adequate for this simple model. 
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A.2.1 The Intake 

The intake compresses the air and reduces its velocity, thereby raising its static temperature. 
In general, the Mach number should be in the range of 0.3 to 0.5 at the compressor inlet. The 
intake-diffusion is assumed to be adiabatic, since die diffusion takes place very rapidly. Hence 
the stagnation temperature remains constant, that is Tt2=Tto- 

Due to imperfect diffusion, a drop in stagnation pressure will generally occur. At high 
velocities the intake efficiency is best specified by a constant kinetic efficiency 


*n 


KE 



(A2) 


For a well designed intake T]ke ~ 0.94 to 0.97, [13,16]. The intake pressure recovery ratio is 
then specified by 
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During off-design operation a fixed-geometry intake would incur spillage drag as a 
consequence of changes in the mass flow, in addition to a reduction in efficiency. However, 
here it is reasonable to assume that for such an advanced vehicle as an aerospace plane, a 
sophisticated variable- geometry intake will be installed, to avoid spillage drag while providing 
the mass flow to match the engine. In other words, in the current model it is assumed that the 
intake itself does not exhibit any off-design behavior (easier stipulated than realized). 


A.2.2 The Compressor 


The compressor provides an additional rise in pressure. The compressor pressure ratio is a 
design parameter. 
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according to [16]. 

Assuming a constant engine speed, there is a constant volume flow through the engine due 
to the constant axial velocity through the compressor [13]. Employing a constant isentropic 
compressor efficiency (r| c =0.96), it can be shown that the off-design pressure ratio is given by 
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where 
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f y > 

f y > 

D = 

y - 1 

e c - 1 

2r lc + e c Y ~ 1 - 1 
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A constant isentropic efficiency has been used as opposed to the more suitable (i.e. more 
realistically constant) polytropic efficiency, in order to be able to obtain a closed-form 
expression for the pressure ratio. 

The compressor exit temperature (Ttf) is determined by 
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The constant volume flow, combined with a constant intake delivery speed of 150 m/s to 
170 m/s [18], allows one to determine the specific mass flow per unit compressor area 



(A7) 


Thus, as will become apparent later, the compressor inlet area constitutes the engine sizing 
parameter which scales the thrust linearly. 

A.2.3 The Combustion Chamber 


In view of its purpose, this model is not concerned with the intricacies of the combustion 
process itself. Instead, simple heating relations for a perfect gas are used with a given heating 
value of the fuel. In this particular case, hydrogen is used in all three propulsion modes. 

For maximum thrust, heat must be added up to the limit set by the TET, approximately 
1600 to 1700 K. For the specific fuel flow one obtains 
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The burning efficiency rib is a correction factor for imperfect combustion, a typical value is 
0.96 [13,17]. H is the heating value of the fuel, H= 144.6 10^ J/kg. An inevitable 
by-product of the heat addition is a loss in stagnation pressure. A reasonable estimate is 3% to 
7%. [13]. 


A.2.4 The Turbine 


In the turbine the gas provides the work needed for compression. 
C P air ( T ,3- T ,2)=’>m‘ : P 8as ( 1+f b)( TCT - T ,5) 
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where rjm represents a mechanical efficiency factor for the turbo-machinery, ri m = 0.99. 
Employing constant isentropic turbine efficiency % the pressure ratio is determined by 
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A.2.5 The Afterburner 

The afterburner increases the specific thrust considerably, an important consideration for a 
boost vehicle like an aerospace plane. Since the afterburner is not followed by any rotating 
machinery the exit temperature (Ttexii) can be higher than the TET. A reasonable value is 2000 
K, [16]. 

The specific fuel flow is given by 
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exit 


(All) 


Again, a slight drop in stagnation pressure is inevitable. 
A.2.6 The Nozzle 


A straightforward sensitivity analysis will show that, for optimal performance, one needs 
to install a convergent-divergent nozzle (also known as a condi nozzle) on a high performance 
vehicle, in order to let the gas expand to near ambient pressure. The decrease in performance, 
using only a convergent nozzle becomes especially noticeable at supersonic speeds. However, 
an upperlimit should, realistically, be imposed on the expansion ratio ER=A ex it/Ath r0 at 1° 
accordance with engine size limitations. 

From the known mass flow through the engine one can determine the necessary throat size 


y + l 



Both the ideal expansion and the given mass flow imply a continuously adjustable nozzle 
geometry. The exit Mach number is determined from either 
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for ideal expansion p e xit = P0> or 
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for maximum expansion, where the maximum expansion ratio (ER) is determined from 


ER = 


exit 


throat 


where A ex it/A c is determined by the engine size and is taken to be approximately 4.0. 
Subsequently, the exit velocity is determined from 
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(A 15) 


where the nozzle velocity factor C v =0.96 is introduced to adjust for nozzle losses. 
The thrust per unit of compressor inlet area is determined from 


= A^exit (‘ + f b +f ah) - V ol + ER ^f~(P«il " P 0 ) 


(A 16) 


and the specific impulse is given by 
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(A 17) 


A.3 THE RAMJET 

The ramjet is more or less a skeleton turbojet with the turbo-machinery taken out since there 
is no need for a compressor, the ram compression by itself is effective enough. The ramjet 
consists of an intake, a combustion chamber and a nozzle. All pertinent equations for these 
components are stated in the previous section on the afterburning turbojet A number of (small) 
changes in some parameter values have been applied, consistent with the different operating 
range of the engine. 

The intake delivery Mach number has been set to a constant value of 0.5. The massflow 
equation (A7) must now be restated in terms of constant intake delivery Mach number 
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The exit temperature has been raised to 2500 K, while the nozzle exit to combustor inlet area 
ratio Aexit/Ac has been raised to 6.0 to boost performance. 


A.4 THE SCRAMJET 


The scramjet is very similar to the ramjet. The main difference is that the combustion takes 
place while the air is moving at supersonic velocities. The reason to apply supersonic 
combustion is to avoid the excessive losses in stagnation pressure, which are introduced first, 
when the air is decelerated by means of an intake to subsonic velocities for combustion, and 
afterwards, when it is accelerated again by means of a nozzle to supersonic velocities. A 
schematic diagram of a scramjet is presented in Figure A2. 


A.4.1 The Intake 


For operation at very high velocities reference [14] suggests using a different type of 
process efficiency, which more nearly stays constant: Kp = 0.93. It can be shown that the 
kinetic efficiency is related to Kd by 


^ KE = K D 0 ~ K d) 
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(A 19) 


and T|ke determines the recovery pressure ratio in the usual way. Additionally, the intake 
delivery Mach numbers are now supersonic. The recommended diffusion ratio M 0 /M 3 is about 
3.0 for optimal fuel specific impulse. 

A.4.2 The Combustion Chamber 

Combustion in a combustor with constant cross section area leads to adequate performance 
at high speeds. However, at the lower speed range constant cross section area combustion 
leads to thermal choking and rapid degradation of performance. To avoid this one needs to 
employ a variable cross section area combustor. A straightforward implementation of such a 
scheme however, would increase the complexity of the model beyond the sought after relative 
simplicity. A particularly elegant way to circumvent this problem is to employ the following 
relation between static pressure and combustor area 

E 

l- E 


(A20) 

Setting E equal to 1.0, one obtains the constant cross section area flow, while setting E equal 
to 0.0 results in constant pressure flow, which in principle cannot choke. By using die above 
relation with linearly increasing E from 0.0 to 1.0 over the Mach number range from 3.0 to 
1 1 . 0 , one can alleviate the performance degradation. 




The equations for one-dimensional flow with heat addition employing (A20) can be found 
in [13]. Given the initial state (M 3 , p t 3 , T t 3 ) and a certain heat input as determined by the 
fuel-to-air ratio fb, the final state (M 5 , p^, T^) can be found. 
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The Mach number at the burner exit can be found from 


(A21) 




2 

E +yM 3 
E +yMj 


2 

1 + 


1 + 


y - 1 


M 


y - 1 


M 


2 

5 

2 

3 


and the stagnation pressure from 
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An optimal fuel-to-air ratio as a function of free-stream Mach number is obtained from 
Swithenbank [Fig. 10, Ref. 14], Therefore the final state of the gas can be determined. In 
case the fuel mixture scheme leads to thermal choking the amount of fuel is reduced to the 
critical amount. The critical fuel-to-air ratio is determined employing (A20) and (A21) in 
reverse order, substituting M 5 = 1.0 into (A21) to obtain T^ and using (A20) next to evaluate 
fb- The supersonic combustion process, i.e. adding heat at supersonic velocities, leads 
unavoidably to a large stagnation pressure loss. However the loss is significantly smaller than 
if subsonic combustion had been used. 


A.4.3 The Nozzle 


In this case, since the flow is already supersonic, the nozzle consists of a divergent section 
only. Again, as with the turbojet and ramjet, the expansion ratio is limited. A value of 6.0 for 
the Aexit/Ac ratio seems reasonable (cf. [15]). The exit Mach number is determined from 
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and the exit velocity is determined from the pertinent equation as stated in the turbojet nozzle 
section. The same remark applies to the specific thrust (here, referenced to the combustor cross 
section area) and the I sp . 



A.5 CONCLUSION 


Again, it should be stressed that the described model is one of relative simplicity and it does 
not occupy itself with the particular details of engine design. However, it adequately describes 
the tendencies of a true multi-mode engine and provides us with a sufficient number of 
parameters to tinker with and refine the modeled engine behavior over a wide range of vehicle 
operating conditions. 
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Figure A1 Basic Components of a Turbojet 
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Figure A2 Basic Components of a Scramjet 




